Corresponding author: 

Dr. Hongjian Feng 

School of Science, Tianjin Polytechnic University, Tianjin 300160, China 

Department of Physics and Astronomy, University of Missouri, Columbia, MO 

65211, USA 

Email address: 

fenghongjian@gmail.com 

fenghongjianl978@gmail.com 



1 



The role of Coulomb and exchange interaction 
on the Dzyaloshinskii-Moriya 
interaction (DMI) in BiFeC>3 



Hongjian Feng a,b 



a School of Science, Tianjin Polytechnic University, Tianjin 300160, China 
b Department of Physics and Astronomy, University of Missouri, Columbia, MO 



65211, USA 



Abstract 

Ab initio calculations show that the Dzyaloshinskii-Moriya interaction(DMI)and net 
magnetization per unit cell in BiFeOs are reduced when U is increasing from to 
2.9 eV, and independent of J. Interestingly, the DMI is even destroyed as U exceeds 
a critical value of 2.9 eV. We propose a simple model to explain this phenomenon 
and present the nature of the rotation of the magnetization corresponding to altered 
antiferrodistortive distortions under DMI in BiFe03. 
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1 Introduction 



Multiferroic materials have attracted much interest due to the coexistence 
of magnetic and ferroelectric ordering in single phase. The coupling of the two 
ordering leads to the so-called magnetoelectric effect in which the magnetiza- 
tion can be tuned by the external electric field, and vice versa [TJ2EEE]- These 
materials have potential applications in information storage, the emerging field 
of spintronics,and sensors. BiFeC>3 is the rare one in nature, which possess both 
weak ferromagnetism and ferroelectric characteristics in single phase [BTrfSfS] . 
It has long been known to be ferroelectric with a Curie temperature of about 
1103 K and antiferromagnetic(AFM) with a Neel temperature of 643 K. The 
Fe magnetic moments are coupled ferromagnetically in (1 1 1) plane and anti- 
ferromagnetically in the adjacent plane along [111] direction, which is known 
as the G-type AFM order. The rhombohedral distorted perovskite structure 
with space group R3c permits a canting of AFM sublattice caused by the 
antisymmetric Dzyaloshinskii-Moriya interaction(DMI), resulting in a weak 
ferromagnetism. However, there is a spiral spin structure in which the AFM 
axis rotates through the crystal with a long-wavelength period of 620A. The 
cancellation of magnetization should be suppressed partly in thin film [B], or 
by partly substitution of magnetic transitional metal ions in B sites, as shown 
in our previous report [10]. It is known that the ferroelectricity in BiFeOa is 
produced by the lone Bi-6s stereochemically active pair induced by the mixing 
between the (ns) 2 ground state and a low- lying (ns) 1 (np) 1 excited state, which 
can only occur if the cation ionic site does not have inversion symmetry, while 
the weak ferromagnetism is mainly attributed to Fe 3+ ions. Therefore the cou- 
pling between the electric and the magnetic ordering becomes weak in BiFe03, 
which agrees with the fact of large difference between the Curie temperature 
and AFM Neel temperature. There exists another structural distortion, so- 
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called antiferrodistortive(AFD) distortion, which is formed by the alternating 
sense of rotation of the oxygen octahedra along [1 1 1] direct ion [TT]. In our pre- 
vious paper, we have shown that the rotation of the oxygen octahedra couples 
with the weak ferromagnetism due to the DMI, using Ab initio calculations 
with considering the spin-orbital(SO) coupling effect and the noncollinear spin 
configuration [12]. In strongly correlated materials, e.g. multiferroics, the on- 
site Coulomb(L r ) and exchange interaction(J) has been proposed to properly 
describe the partly filled localized d orbitals within density functional the- 
ory (DFT). One may wonder whether U and J will have an impact on the 
magnetization through DMI taking into account SO interaction. How do these 
parameters influence the DMI, and further the magnetization? What is the 
origin of coupling between the rotation of oxygen octahedra and the resulting 
magnetization in terms of DMI in BiFeOs? In this paper we have proposed 
a transparent physical interpretation for the abovementioned questions,using 
first-principles calculations based on the DFT. 

The remainder of this article is structured as followsdn section 2, we 
presented the computational details of our calculations. We provided the cal- 
culated results and discussions in section 3. In section4, the conclusion based 
on our calculation were given. 

2 Computational details 

Our calculations were performed within the local spin density approx- 
imation(LSDA) to DFT using the ABINIT package [13114] . The ion-electron 
interaction was modeled by the projector augmented wave (PAW) potentials 
[T5][T6] with a uniform energy cutoff of 500 eV. Bi 5d, 6s, and 6p electrons, Fe 
4s, 4p,and 3d electrons, and O 2s and 2p electrons were considered as valence 



4 



states. Two partial waves per / quantum number were used. The cutoff radii for 
the partial waves for Bi, Fe, and were 2.5, 2.3, 1.1 a.u., respectively. 6x6x6 
Monkhorst-Pack sampling of the Brillouin zone were used for all calculations. 
We calculated the net magnetization per unit cell and the electronic proper- 
ties within the LSDA+U method where the strong coulomb repulsion between 
localized d states has been considered by adding a Hubbard-like term to the 
effective potential [T7|18f 19] , The effective Hubbard parameter, the difference 
between the Hubbard parameter U and the exchange interaction J (U — J), 
was changing in the range between and 6 eV for the Fe d states. For the same 
value of (U—J), J was varying as 0,0.5, 0.8, and 1 eV, respectively. Taking into 
account the SO interaction, we introduced the noncollinear spin configuration 
to construct the G-type AFM magnetic order with the AFM axis being along 
the x axis in Cartesian coordinates in our Ab initio calculation. 

3 Results and discussion 

In ref. 20 the author suggest that the inversion centers between adjacent 
B sites in ABO3 perovskite structure are destroyed by the displacement of 
the oxygen anions located at the midpoints between them, while the space in- 
version centers between A sites still remains. Therefore ABO3 structure with 
magnetic ions in A sites, such as FeTiOs, should possess a strong coupling be- 
tween the ferroelectric distortions and magnetization. It can not be achieved 
in ABO3 structure with magnetic ions in B sites, such as BiFeOs. That is to 
say the coupling between the ferroelectric distortions and magnetization in 
it shall be neglected. However in BiFeOa there exists another kind displace- 
ment, known as antiferrodistortive(AFD) distortions, caused by the rotation 
of the neighboring oxygen octahedra. Through antisymmetric superexchange 
interaction this AFD displacement couples weakly to the magnetization. In 
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this paper we mainly concentrate on the coupling associated with the DMI 
between the AFD distortions and the magnetization per unit cell. 



For the AFD motion, a rotational vector R has been introduced to de- 
scribe the direction of the rotation of the oxygen octahedra|12j. From Fig. 1, 
the anticlockwise rotation of upper oxygen octahedra and clockwise rotation of 
lower oxygen octahedra correspond to the outward state defined as Ti out . The 
opposite state is defined as Rj n . The rotational angle is 10°in the Cartesian 
coordinates 1 121. 



In our LSDA+U calculation, U and J are defined as 
1 



U 



— - J2 < m,m'\V e e\m,m' >= F°, (1) 



(21 , 



1 ^ , lTr , , . F 2 + F A 

2/(2/ 



J— — — < m, m'\V ee \m, m! >= , (2) 

1) _f-t , 14 ' v ; 



where V ee are the screened Coulomb interaction among the nl electrons. F°, 
F 2 , and-F 4 are the radial Slater integrals for d electrons in Fe. 

The net magnetization per unit cell with respect to Rj n and R OU i in 
Cartesian coordinates for different U and J were listed in table 1. It can be 
seen that J value have no effect on the resulting magnetization when U remains 
constant. For the sake of clarity, only the results obtained with different J value 
for U=0 and 2.9 eV were given in the table . The AFM vector in Cartesian 
coordinates with varying effective Hubbard U were illustrated in Fig. 2, where 
[1 1 1] direction is taken as the z axis as shown schematically in Fig. 3, and 
the x, y,and z component of the magnetization is denoted by M x , M y , and 
M z in the coordinates, respectively. Fig. 3 shows the coupling between the 
rotation of oxygen octahedra and the resulting magnetization per unit cell. 
The arrow indicate the spin direction of Fe for different states. The upper 
section corresponds to the Rj n rotational state, and lower section, the R ut 
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rotational state. M to t a i is the net magnetization per unit cell. The dashed line 
arrow is the unstable rotational state. [1 1 1] crystal direction is selected as the 
z axis, and the AFM order is arranged along x axis. It is clearly shown that the 
easy axis of the magnetization is y axis when G-type AFM order is arranged 
along the x axis, taking into account the SO interaction and the unconstrained 
freedom of spin. The antisymmetric interaction of the neighboring Fel and 
Fe2 ions leads to the canting of the magnetic moment of them away from 
their original direction(x axis) and a resulting magnetization mainly in y axis, 
which arises from the DMI only occurring when the inversion symmetry is 
broken. As U is approaching from eV to 2.9 eV,the net magnetization is 
reversed by the opposite rotation of the oxygen octahedra in terms of the 
reversal of M y , and decreases with increasing of U. However, M y does not 
change sign with the altered AFD motion when U exceeds a critical value of 
2.9 eV, say 3 eV, implying that the net magnetization only deviates slightly 
from the original direction and does not experience a significant rotational 
angle greater than 90°. As U attain to be the critical value, the DMI caused 
by the antisymmetric superexchange interaction is eliminated with the strong 
on-site Coulomb interaction. The AFD distortions do not couple with the 
magnetization. 

In order to obtain an unambiguous interpretation for the effect of Coulomb 
and exchange interaction on the net magnetization, we need to recap the 
DMI on the coupling of neighboring Fel and Fe2 sites. We have for the in- 
teraction of neighboring Fel and Fe2 sites by the second order perturbation 
calculation [2TI22I23] 



E fIi,Fc2 - J?el,Fe2(Sl ' S 2 ) + D^el^^l X S 2 ) + S(R) ■ r^ 1Fe2 ■ S 2 . (3) 
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The first term on the right hand side of the Eq. (3) corresponds to the usual 
isotropic superexchange interaction, and the second term is the DMI. Provided 
the long range pseudodipolar interaction is neglected, we get the Hamiltonian 
for the system 

HsiFeOa = ~ 2 Jlt,2j"Slj ' $>2j + ^ Dli,2j'Sij X (4) 

The first term comes from the symmetric superexchange, and the second one 
is the antisymmetric DMI contribution. Ju^j in the first term is a constant 
similar to the exchange interaction, and does not contribute to the DMI. This 
can account well for our calculated results that the exchange parameter J has 
nearly no effect on the rotation of the magnetization. D is the DMI constant 
associated with the crystal field and determined by the sense of rotation of 
the neighboring oxygen octahedra(Rj n or R OU i). D reads by the second order 
perturbation in the case of one electron per ion 

DF 2 ei,Fe2 = (4:i/U)[b nn >(Fel—Fe2)C n r n (Fe2—Fel)—C nn '(Fel—Fe2)b n ' n (Fe2—Fel)], 

(5) 

where U is the energy required to transfer one electron from one site to its 
nearest neighbor, a parameter similar to on-site Coulomb interaction in our Ab 
initio computation, and inversely proportional to D. This is consistent with 
our calculated results that the absolute value of net magnetization is inversely 
proportional to the Hubbard parameter U. Magnetization does not reverse its 
direction in terms of the changing of the AFD displacement, especially when 
U is greater than the critical value of 2.9 eV , this indicates that in this case U 
is large enough to make DMI being disappeared. We have also calculated the 
band gap for different U corresponding to H in and R OM t, respectively. From 
Fig. 4, it can be seen that the curve becomes relatively flat when U reaches 
the critical value of 2.9 eV. Thereafter, we chose this value to describe the 
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electronic property in the following. It is worth mentioning that the band gap 
to R in is greater than to R ou t, indicating that the AFD motion corresponding 
to R OU i tend to reduce the crystal-field splitting , and consequently the band 
gap. 

In order to analyze the rotation of magnetization under DMI, we have 
calculated the Orbital-resolved density of states(ODOS) for Fel and Fe2 cor- 
responding to R iTl and H ut m Fig. 5, Fig. 6, Fig. 7, and Fig 8, respectively. 
Fig. 5 is the ODOS for Fel corresponding to Rj n rotational state. The vertical 
line indicates the Fermi level. All the states occupied in the valence band are 
spin- up electrons (majority spin as defined). It means the spin direction for Fel 
to R in is positive. Fig. 6 is the ODOS for Fe2 corresponding to H in rotational 
state. The vertical line indicates the Fermi level. All the states occupied in the 
valence band are spin-down electrons(majority spin as defined). It means the 
spin direction for Fe2 to R« n is negative. Fig. 7 is the ODOS for Fel corre- 
sponding to H ut rotational state. The vertical line indicates the Fermi level. 
All the states occupied in the valence band are spin-up electrons (majority 
spin as defined). It means the spin direction for Fel to R OU i is positive. Fig. 
8 is the ODOS for Fe2 corresponding to H ou t rotational state. The vertical 
line indicates the Fermi level. All the states occupied in the valence band are 
spin-down electrons (majority spin as defined). It means the spin direction for 
Fe2 to Hout is negative. Let us come back to Fig. 3. In order to make the net 
magnetization reversed, the magnetic moment of Fel and Fe2 can rotate from 
the original direction corresponding to R« n (Fig. 3(a)) either to the dashed 
line arrow(Fig.3 (b)) required greater energy barrier, or to the real line ar- 
row required smaller energy barrier. From Fig. 5 to Fig. 8, one can see that 
the spin-up electrons in the occupied valence band for Fel and the spin-down 
electrons in the occupied valence band for Fe2 do not change their in-built spin 
direction when rotational vector is changing from R» n to H ut- This confirms 
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that the spin direction of Fel and Fe2 only deviates slightly from the initial 
states to the final states as shown in the real line arrow in Fig. 3 (b) corre- 
sponding to R ut. It is worth pointing out that d x 2_ y 2 orbital for Fel and Fe2 
is split from the doubly degenerate e g states and tend to overlap with d xy , d yz , 
and d xz orbitals in the triply degenerate ti g states, indicating that the AFM 
DMI is mainly attributed to the e g -e g AFM interaction which is greater than 
the t2g-t2 9 AFM interaction. 

4 Conclusion 

Magnetization can be reversed by the altering sense of rotation of the 
oxygen octahedra in BiFeC>3 when U is smaller than the critical value of 2.9 
eV, and the absolute value of magnetization is decreasing as U is ranging from 
to 2.9 eV. Magnetization does not reverse with altered AFD displacement 
when U exceeds the critical value, indicating that the DMI is even prohibited 
in this case. The rotation of magnetization is fulfilled by slight deviation of 
magnetic moment of Fel and Fe2 around rather than reversal of them. 
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Table 1 

Magnetization per unit cell with respect to different value of U andJ. 

U(eV) 0.5 0.8 1 1 

J(eV) 0.5 0.8 1 

E^in I*out -f^in Rout ^Mii R-out -f^in Rout -f^in Rout 

M x (hb) 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 

M y (fj, B ) 0.4259 -0.0812 0.4259 -0.0812 0.4259 -0.0812 0.4259 -0.0812 0.0351 -0.0679 

M z (fM B ) -0.1013 0.0000 -0.1013 0.0000 -0.1013 0.0000 -0.1013 0.0000 -0.0493 0.0056 

U(eV) 2 2.5 2.8 2.9 3.4 

J{eV) 0.5 

Rin Rout Rin Rout Rin Rout Rin Rout Rin Rout 

M x (fj, B ) 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 

M y (fj, B ) 0.0416 -0.0337 0.0325 -0.0251 0.0365 -0.0188 0.0313 -0.0176 0.0313 -0.0176 

M z (hb) -0.0408 0.0108 -0.0366 0.0147 -0.0431 0.017 -0.0406 0.0168 -0.0406 0.0168 

U(eV) 3.7 3 4 5 6 

J{eV) 0.8 

Rin Rout Rin Rout Rin Rout Rin Rout Rin Rout 

M x (fi B ) 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 

M y {n B ) 0.0313 -0.0176 0.0237 0.0012 0.0172 0.0033 0.0176 0.0085 0.0178 0.0111 

M z (fM B ) -0.0406 0.0168 -0.0283 -0.0049 -0.0249 -0.0055 -0.0186 -0.0053 -0.0157 -0.0025 
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Figure captions: 

Fig.l The rotational vectors of the AFD distortions in BiFeC>3. The shaded 
cage denote the oxygen octahedra, and Fe is inside the cage. 

Fig. 2 AFM vectors with respect to U. 

Fig. 3 Schematic diagram for the coupling between the rotation of oxygen 
octahedra and the resulting magnetization in unit cell in BiFeC>3. The arrow 
denote the direction of magnetization. 

Fig. 4 Band gap for R in and R OU i with respect to U. 

Fig. 5 ODOS for Fel d xy ,d yz ,d z 2,d xz , and d x 2_ y 2 orbitals to R in . 

Fig. 6 ODOS for Fe2 d xy ,d yz ,d z 2,d xz , and d x 2_ y 2 orbitals to R in . 

Fig. 7 ODOS for Fel d xy ,d yz ,d z 2,d xz , and d x 2_ y 2 orbitals to R ou t. 

Fig. 8 ODOS for Fe2 d xy ,d yz ,d z 2,d xz , and d x 2_ y 2 orbitals to H ou t- 
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Fig. 2. AFM vectors with respect to U. 
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Fig. 3. Schematic diagram for the coupling between the rotation of oxygen octahe- 
dra and the resulting magnetization in unit cell in BiFe03. The arrow denote the 
direction of magnetization. 



16 




Fig. 4. Band gap for Rj n and R ou t with respect to U. 
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Fig. 5. ODOS for Fel d xy ,d yz ,d z 2,d xz , and d x 2_ y 2 orbitals to 
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E-E F (eV) 

Fig. 6. ODOS for Fe2 d xy ,d yz ,d z 2,d xz , and d x 2_ y 2 orbitals to 
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Fig. 7. ODOS for Fel d xy ,d yz ,d z 2,d xz , and d x 2_ y 2 orbitals to R, 



out- 
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Fig. 8. ODOS for Fe2 d xy ,d yz ,d z 2,d xz , and d x 2_ y 2 orbitals to Rout- 
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